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ABSTRACT 

The  nonlinear  instability  of  the  boundary  layer  on  a  heated  flat  plate  placed  in  an  on¬ 
coming  flow  is  investigated.  Such  flows  are  unstable  to  stationary  vortex  instabilities  and 
inviscid  traveling  wave  disturbances  governed  by  the  Taylor-Goldstein  equation.  For  small 
temperature  differences  the  Taylor-Goldstein  equation  reduces  to  Rayleigh’s  equation.  When 
the  temperature  difference  between  the  wall  and  free  stream  is  small  the  preferred  mode  of 
instability  is  a  streamwise  vortex.  It  is  shown  in  this  case  that  the  vortex,  assumed  to  be  of 
small  wavelength,  restructures  the  underlying  mean  flow  to  produce  a  profile  which  can  be 
massively  unstable  to  inviscid  traveling  waves.  The  mean  state  is  shown  to  be  destabilized 
or  stabilized  to  inviscid  waves  depending  on  whether  the  Prandtl  number  is  less  or  greater 
than  unity. 
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1.  Introduction 


Our  concern  is  with  the  nonlinear  instability  of  forced  convection  boundary  layers  over 
horizontal  heated  flat  plates  and  the  subsequent  secondary  instabilities  of  these  flows.  In 
a  previous  paper,  Hall  and  Morris  (1992),  hereafter  referred  to  as  HM,  considered  the  lin¬ 
ear  aspects  of  longitudinal  vortex  instabilities.  In  general  it  was  shown  that  this  convec¬ 
tive  instability  develops  in  a  nonparallel  manner  and  cannot  be  adequately  described  by  a 
quasi-parallel  stability  theory  of  the  type  discussed  by  for  example  Wu  and  Cheng  (1976), 
Moutsoglou,  Chen  and  Cheng  (1984),  Chen  and  Cheng  (1981).  At  higher  values  of  the  con¬ 
trolling  stability  parameter,  the  Grashof  number,  nonparallel  effects  can  be  taken  care  of  in 
a  self-consistent  asymptotic  manner  provided  that  the  vortex  wavelength  is  not  too  large. 
In  the  present  paper  we  will  be  concerned  with  extending  the  work  of  HM  into  the  strongly 
nonlinear  regime,  in  addition  we  will  be  concerned  with  the  subsequent  three-dimensional 
unsteady  breakdown  of  the  flow  induced  by  large  amplitude  vortex  structures. 

In  addition  to  longitudinal  vortex  structures  both  Tollmien-Schlichting  and  inviscid  waves 
are  possible  causes  of  instability  in  a  heated  boundary  layer.  If  buoyancy  forces  are  not  too 
large,  then  the  inviscid  inodes  are  found  to  satisfy  Rayleigh’s  equation,  otherwise  they  satisfy 
the  Taylor-Goldstein  equation.  We  show  that,  when  buoyancy  forces  are  sufficiently  large 
enough  to  alter  the  zeroth  order  inviscid  instability  problem,  they  also  enter  the  equations 
for  the  basic  state  and  cause  the  mean  velocity  and  temperature  fields  to  be  coupled. 

In  order  to  make  progress  with  an  analytical  solution  of  the  strongly  nonlinear  instability 
problem  associated  with  longitudinal  vortex  instabilities,  we  use  the  asymptotic  structure 
given  by  Hall  and  Lakin  (1988)  in  the  context  of  small  wavelength  Gortler  vortices.  In  this 
limit,  which  corresponds  to  the  far  downstream  behavior  of  a  fixed  wavelength  vortex,  the 
vortices  are  confined  to  a  finite  part  of  the  boundary  layer  and  indeed  where  they  exist  the 
mean  state  adjusts  to  allow  them  to  remain  neutral.  We  show  that  in  forced  convectivion 
boundary  layers  the  vortices  occupy  a  region  adjacent  to  the  wall  and  that  the  thickness  of 
this  region  increases  linearly  with  Grashof  number.  The  results  we  present  are  for  similarity 
solutions  of  the  nonlinear  vortex-mean  flow  interaction  problem  but  previous  experience  with 
the  related  Gortler  problem,  Hall  (1988),  Hall  and  Lakin  (1988)  suggests  that  the  structure 
of  the  nonsimilar  solutions  is  not  significantly  different. 

A  significant  result  which  we  find  is  that  for  Prandtl  numbers  less  than  unity  favorable 
pressure  gradient  flows  can  have  inflectional  streamwise  mean  velocity  profiles  induced  by 
vortices.  This  suggests  that  a  major  role  of  vortex  instabilities  might  be  to  modify  favorable 
pressure  gradient  flows  to  make  them  unstable  to  the  relatively  much  more  dangerous  inviscid 
instabilities.  We  investigate  the  Rayleigh  instability  problem  for  such  flows  and  show  that 
vortices  induce  a  class  of  unstable  Rayleigh  waves  over  a  large  range  of  wavenumbers. 
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The  procedure  adopted  in  the  rest  of  this  paper  is  as  follows:  in  §2  we  formulate  the 
nondimensional  equations  governing  vortex  instabilities  in  heated  boundary  layers.  In  §3 
these  equations  are  solved  for  small  vortex  wavelengths  and  in  §4  a  similarity  solution  is 
obtained.  Some  large  Grashof  number  properties  of  the  similarity  solution  are  discussed  in 
§5  where  we  show  the  importance  of  the  Prandtl  number  in  determining  the  shape  of  the 
induced  mean  flow.  In  §6  we  discuss  the  inviscid  instability  of  the  vortex  states  discussed  in 
§3,  4,  5  and  finally  in  §7  we  draw  some  conclusions. 

2.  Formulation 

Consider  the  viscous  flow  over  a  semi-infinite  flat  plate  with  a  typical  lengthscale  L  in 
the  flow  direction.  Suppose  that  the  fluid  speed  at  infinity  is  UoUe(x* / L)  where  x*  measures 
distance  along  the  wall.  We  define  the  parameters  R  and  G  by 

R  =  U0Lu~\  G  =  L3g(T0  -  T00)/3i^~2R~:i/2.  (2.1a,  b) 

Here  v  is  the  kinematic  viscosity,  g  the  acceleration  due  to  gravity  and  /?  is  the  coefficient  of 
expansion  of  the  fluid.  The  temperatures  T0  and  respectively  denote  a  constant  reference 
temperature  and  the  fluid  temperature  in  the  freestream.  We  concern  ourselves  with  the  limit 
R  — »  oo  with  G  held  fixed  and  determine  the  strongly  nonlinear  vortex  flows  which  occur  in 
this  limit. 

Suppose  next  that  the  wall  temperature  is  given  by 

T  =  Too  +  {Too  —  T0)T(xm/L)  (2.2) 

and  that  we  define  dimensionless  variables  (x,y,z)  by 

(x,y,z)  =  (x-L-\y'L-lRl'\z’L-lRl/2),  (2.3) 

where  ym,zm  denote  distance  in  the  normal  and  spanwise  directions  respectively.  We  define 
a  dimensionless  velocity  field  ( u,v,w )  associated  with  ( x,y,z )  by  writing 

(u*,  e*,  «;*)  =  U0(u,  R~l'2v,  R~1/2w).  (2.4) 

The  pressure  pm  is  then  written  in  the  form 

P“  =  PL  o  jpoU)  +  PiU)/#172  +  +  .  ..j  ,  (2.5) 

where  p  is  a  typical  fluid  density.  The  fluid  density  is  given  by 

p-  =  pli  +  Mr  -  rj]  (2.6) 
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and  we  shall  make  the  Boussinesq  approximation  throughout  this  work  so  that  changes  in 
p  can  be  ignored  unless  they  are  multiplied  by  gravity.  The  temperature  field  is  then  made 
dimensionless  by  writing 

T’  =  Too  +  (To  -  Too)T.  (2.7) 

In  the  Boussinesq  approximation  the  continuity  equation,  momentum  equations  and  energy 
equations  are: 

V  •  u  =  0, 

\ 

^•V)^  =  (-pyG  +  GT)-py  +  (2.8a,  b,c) 

~Vz 


(m-V)T  =  -V2T, 

<7 

where  a  is  the  Prandtl  number,  V  =  (dx,dy,dx)  and  V2  =  (dy  +  d]  +  R~ldl).  The  boundary 
conditions  appropriate  to  (2.8)  are 


u  =  0,  T  =  T(x),  y  -  0, 
u  — »  ue?  — »  0,  T  — *  0,  y  — ►  oo. 


(2.9a,  b) 


Note  here  that  if  the  instability  is  caused  by  for  example  nonuniform  wall  heating  the  con¬ 
ditions  at  the  wall  should  be  modified  as  in  HM.  In  the  absence  of  any  instabilities  we  write 
u  =  (u,v,  0),  T  =  T  and  the  basic  state  is  determined  in  the  limit  R  —*  oo  by 


ux  +  vy  =  0, 

uux  +  vuy  =  ueuex  +  uyy  -  Spx 
Py  =  T , 

uTx  +  vTy  =  I”. 

a 

u  =  v  =  0,  T  =  T,y  =  0, 
u  — >  ue,  T  — ►  0,  y  — >  oc. 


(2.10a -f) 


Here  we  have  set  p0x  =  —ueuex  and  introduced  the  buoyancy  parameter  S  =  G/R  which  we 
take  to  be  0(1).  The  longitudinal  vortex  structure  of  HM  first  occurs  for  G  =  0 (R°)  but  our 


ultimate  concern  is  with  the  interaction  of  vortex  structures  and  Taylor- Goldstein  inviscid 
waves  which  exist  for  G/R  ~  0 (R°).  Hence  it  is  convenient  at  this  stage  to  allow  for  a 
‘buoyancy’  coupling  between  the  x,y  momentum  equations,  even  though  most  of  the  results 
given  in  this  paper  will  be  for  the  case  5  =  0.  Note  also  that  (2.10c)  can  be  integrated  and 
the  result  substituted  into  (2.10b)  to  give 

uux  +  vUy  —  ueuex  +  uvy  +  S  f  Txdy.  (2.11) 

J  OO 

This  coupling  has  a  significant  effect  on  similarity  solutions  for  the  basic  state.  In  particular 
if  we  choose  ue  =  xn  then  the  wall  temperature  T  must  be  proportional  to  i5"/2-1/2  for  a 
similarity  solution  to  exist. 

In  HM  the  linearized  version  of  (2.8)  was  solved  for  the  case  n  =  1;  it  was  shown  that 
at  0(1)  values  of  G  the  longitudinal  vortex  instability  induced  by  heating  develops  in  a 
nonparallel  manner.  At  higher  values  of  G  interest  centres  on  the  right-hand  branch  of  the 
neutral  curve  where  vortices  of  wavenumber  a  are  neutral  for  G  ~  a4.  In  fact  this  high 
wavenumber  regime  is  important  because  a  disturbance  of  fixed  physical  wavelength  will 
correspond  to  large  values  of  the  local  Grashof  and  wavenumbers  at  sufficiently  large  values 
of  x.  Thus  the  regime  is  particularly  relevant  to  experiments  where  the  vortex  instability  is 
caused  by  very  small  background  disturbances  which  must  evolve  over  a  significant  distance 
before  they  develop  in  a  nonlinear  manner.  It  is  with  the  latter  nonlinear  problem  which  we 
will  concern  ourselves  in  this  paper. 

3.  The  Strongly  Nonlinear  Evolution  Equations  for  Small  Wavelength 
Streamwise  Vortices 

The  nonlinear  structure  which  we  develop  in  this  section  is  based  on  the  related  analysis 
of  Hall  and  Lakin  (1988)  for  Gortler  vortex  flows.  The  major  difference  between  the  problems 
however  is  that  in  the  convection  problem  small  wavelength  neutral  vortices  occur  near  to  the 
boundary  whilst  in  the  Gortler  problem  the  vortices  move  away  from  the  wall  to  the  position 
in  the  flow  where  Rayleigh’s  criterion  is  most  violated.  This  difference  leads  to  significant 
differences  in  the  corresponding  nonlinear  structures.  In  fact  the  heated  boundary  layer 
has  some  similarities  with  the  strongly  nonlinear  Taylor  vortex  problem,  Denier  (1992). 
However,  the  fact  that  the  Taylor  instability  occurs  in  a  bounded  region  causes  even  the 
Taylor  and  convection  problems  to  differ  significantly  at  high  values  of  the  appropriate 
control  parameters. 

In  the  convection  problem,  small  wavelength  vortices  feel  the  local  strength  of  the  desta¬ 
bilizing  bouyancy  force  through  |^,  the  vertical  gradient  of  the  basic  temperature  field.  In 
general,  this  gradient  is  a  maximum  at  the  wall  so  the  instability  is  initiated  there.  In  fact 
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it  was  shown  in  HM  that  in  the  neutral  case  the  dominant  balance  is  between  spanwise 
diffusion  and  buoyancy  effects  in  the  y  momentum  equation  and  between  spanwise  diffusion 
and  the  convective  operator  in  the  energy  equation.  It  is  this  balance  which  leads  us  to  the 
appropriate  generalization  of  Hall  and  Lakin  (1988)  to  the  present  problem. 

Suppose  then  that  a  vortex  of  wavelength  a  exists  in  a  region  of  depth  0(1)  adjacent 
to  y  =  0.  Since  we  have  assumed  that  G/R  is  0(1)  then  we  must  consider  the  case  when 
a  =  0{R1'4).  Thus  we  write 

G  =  Ga4  +  •  •  • 

(3.1a, b) 

with  G/R  =  S  +  •  •  •  . 

In  order  to  recover  the  situation  when  buoyancy  is  not  important  we  must  then  consider  the 
second  limit  S  — »  0.  The  neutral  right  hand  branch  modes  for  a  >  1  have  temperature, 
pressure  and  velocity  perturbations  0,p,u,v,w  such  that 


v  ,  v  ,  v 
~  0(a2),  —  ~  0(1;,  r  ~ 


0(a2),  4 

in 


0(a). 


(3.2) 


6  P  u 

The  first  of  the  above  scalings  follows  from  our  earlier  remarks  concerning  the  balances 
GO  ~  v22 ,  022  ~  vTy  whilst  the  remaining  balances  follow  from  the  choice  of  balances 
Pz  ~  u;„,  u22  ~  vuy,  vy  ~  w2  in  the  z,x  momentum  equations  and  the  continuity  equation 
respectively.  We  then  fix  the  overall  size  of  the  disturbance  so  that  the  mean  flow  and 
temperature  field  corrections  driven  by  the  disturbance  are  comparable  with  the  unperturbed 
flow.  Thus  we  require  that  u  •  Vu  ~  0(1)  ~  ii  ■  V0;  we  therefore  expand  u,v,w,p  and  T 
appearing  in  (2.8)  in  the  form 

u  =  u0  +  y)  H - b  a~l[U0(x,  y)E  +  C.C.]  +  •  •  • 


v  =  v0  +  a  2/%i  (x,j/)  + - b  a[V0(x,y)E  -f  C.C.J  + 

w  =  [W0(x,y)E  +  C.C.  ]  +  ••• 


(3.3a,  b,c,  d,e) 


p  =  [ a2P0(x,y)E  +  C.C.}  +  ••• 

T  =  To(x,y)  +  a-V3Tl(Xly)  +  ---  +  a-l[0oE  +  C.C.}  +  ---  . 

In  the  above  expressions  the  function  2 E  =  e'az ,  C.C.  denotes  complex  conjugate  and  •  •  • 
denotes  terms  smaller  (in  terms  of  a)  than  those  immediately  before  this  symbol.  Notice 
also  that  the  mean  (in  terms  of  z)  part  of  the  pressure  must  be  expanded  as 

P  =  %ix,  y)  +  a'2/39i(x,  y)  +  •  ■  •  (3.3f) 

Furthermore,  we  have  anticipated  that  the  correction  terms  in  (3.3)  to  the  mean  and  funda¬ 
mental  are  0(a_2//3)  smaller  than  the  dominant  terms.  This  choice  can  be  inferred  from  the 
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fact  that  the  correction  terms  to  the  mean  must  be  comparable  to  the  depth  of  the  transition 
layer  in  which  the  vortices  decay  to  zero.  Since  the  depth  of  this  layer  is  governed  by  the 
scalings  of  the  linear  neutral  problem  it  must  be  of  depth  0(a-2/3),  see  Hall  and  Lakin  (1988). 


Equations  for  the  vortex 


If  we  substitute  (3.3)  into  (2.8)  and  retain  the  leading  order  fundamental  terms  we  obtain: 


Voy  +  iW o  —  0, 

U0  +  VoUQy  =  0, 


GOq  —  Vo  =  0, 

Wo  ipQ  =  0, 


(3.4a,  b,c,d,e) 


- f  VoToy  —  0. 

<7 

These  equations  determine  only  £/0,  Vo,  Wo,  and  P0  in  terms  of  0O,  and  in  fact  the  equations 
are  only  consistent  if 

eGT  Qy  +  1=0. 


We  integrate  this  equation  to  give 

To  =  T( i)  -  4-,  (3.41) 

C jG 

where  we  have  satisfied  the  required  boundary  condition  on  T0  at  y  —  0.  The  temperature 
profile  (3.4f)  is  that  temperature  distribution  which  enables  a  vortex  of  wavenumber  a  1 
to  remain  neutral.  Without  any  loss  of  generality  we  will  now  assume  that  the  fundamental 
temperature  disturbance  do  is  real. 


Equations  for  the  mean: 

If  we  substitute  the  expansions  (3.3)  into  (2.8)  and  equate  the  leading  order  mean  terms 
in  the  continuity,  x,  y  momentum  equations  and  the  energy  equations  and  use  (3.3)  we  obtain 

Uqx  +  t'Uy  =  0, 

UoUox  +  UoUoy 

%y 

UcTor  +  I'0 1  Oy 


—  T  +  u0yy  +  n  ^u0j|V0  |  , 


=  To 


(3.5a,  b,  c,  d) 


+  i>}(. 
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The  (/-momentum  equation  above  can  be  integrated  directly  to  give 

%  =  [VT0dy  +  q0(x,Q).  (3.5e) 

Jo 

The  function  To  appearing  in  (3.5c)  can  then  be  replaced  by  1  —  4^  and  then  (3.5a,b,d)  can 
be  integrated  to  find  Vq2,  tZo,Uo-  It  is  instructive  at  this  stage  to  seek  a  similarity  solution  of 
(3.5)  in  order  to  gain  some  insight  into  the  nonlinear  structure  we  have  obtained. 

We  recall  that  for  the  Falkner-Skan  profile  with  ue  ~  xn  buoyancy  forces  can  only  be 
retained  within  the  similarity  solution  structure  if  T  ~  x^F1.  However  (3.4c)  implies  that  a 
similarity  solution  in  the  presence  of  vortices  is  possible  only  if  T  ~  x^1-n^2.  Thus  a  similarity 
solution  in  the  presence  of  both  buoyancy  effects  and  vortices  is  possible  only  if  n  =  1/3.  In 
fact  the  case  n  =  1/3  is  of  particular  interest  in  vortex- wave  interaction  theory,  since,  as  will 
be  remarked  upon  later,  the  only  possible  similarity  solution  when  small  wavelength  vortices 
interact  with  Rayleigh  or  Taylor- Goldstein  waves  has  n  =  1/3.  In  the  next  section  we  will 
discuss  the  n  =  1/3  case  in  detail,  before  doing  so  we  will  complete  the  description  of  the 
flowfield  for  the  more  general  non  self-similar  case.  However  we  stress  that  the  results  which 
we  obtain  for  this  particular  choice  of  n  are  typical  of  the  other  similarity  solutions. 

The  equations  (3.5a,b,d)  must  in  general  be  integrated  numerically  with  qQ  given  by 
(3.5c)  and  T0  given  by  (3.4f).  However,  we  can  think  of  (3.5d)  as  a  first  order  equation  in  y 
for  Tq2  and  we  formally  write  the  solution  in  the  form 

VS  =  K?(*, 0)  -  2 G  r  \u0%  -  -Fl  iy.  (3.6) 

JO  L  Go-J 

The  function  Vo(x,0)  remains  unknown  at  this  stage  and  for  small  values  of  y  the  integrand 
behaves  like  u0y{x,y)yTx  so  that  if  we  assume  there  is  no  reversed  flow  then  this  quantity 
is  positive  if  Tx  >  0.  Thus  for  similarity  solutions  with  n  <  1,  Vq2  is  a  positive  decreasing 
function  near  the  wall  and  we  expect  that,  since  we  do  not  anticipate  the  presence  of  vortices 
everywhere,  the  vortex  will  vanish  at  some  value  of  the  similarity  variable  r/.  In  fact  for 
similarity  solutions  with  values  of  n  not  satisfying  this  inequality  the  integrand  becomes 
positive  for  large  enough  r/,  thus  the  only  major  change  is  that  the  maximum  of  V0  occurs 
away  from  the  wall.  We  also  expect  that  (3.6)  will  determine  Vq  in  a  finite  range  of  valuesfor 
y  in  the  non  self-similar  case. 

The  solution  (3.6)  is  therefore  valid  for  0  <  y  <  y  and  more  precisely  near  y ,  Vo  ~  |y  —  y\ 
and  viscous  effects  come  back  into  play.  In  fact  a  layer  of  depth  a-2^3  is  required  at  y  in 
order  to  allow  Vo  to  adjust  from  algebraic  decay  for  y  —  y  >>  a~2^3  to  exponential  decay  for 
y  —  y  >>  a-2/3.  The  required  structure  in  this  layer  is  virtually  identical  to  that  given  by 
Hall  and  Lakin  (1988)  and  so  will  not  be  repeated  here.  It  suffices  to  say  that  Vo  satisfies 
the  second  Painleve  equation  in  this  layer  and  that  across  the  layer  Uo,u0y,v0,To,Toy  and  qo 


are  continuous.  Thus  above  the  transition  layer  we  retain  the  expansions  (3.3)  but  with  all 
the  2  dependent  terms  set  equal  to  zero  so  that  the  leading  order  problem  for  Uq,vq,  etc.  in 
this  upper  layer  is  identical  to  that  for  the  unperturbed  flow.  We  can  therefore  write  down 
the  following  ‘composite’  problem  for  the  whole  flow  field. 

UO*  +  voy 

UqUOx  +  VQllOy 


UqT  QX  +  VqT  Qy 

where  H  is  the  Heaviside  function 

'1,  s  >  0 

H(s)  =  < 

.0,  s  <  0. 

The  appropriate  boundary  conditions  are 

Vo  =  Vo(x,0),uo  =  Vo  =  0,  Tq  -  T,  y  =  0 

u0  — * >  ue,  To  —> ►  0,  y  — +  oo  (3.8) 

and  q0l(x,  0)  =  —  /0°°  T0xdy. 

We  note  that  at  the  point  y  where  Vq  =  0  the  functions  u0,uoy,vo,To,T0y  are  to  be  made 
continuous  and,  since  (Ho,  To)  — >  (ue,0)  as  y  — *  oo,  then  Vo(x,0),  HOy(x,0),  must  adjust 
themselves  in  order  that  this  limit  is  achieved.  Note  also  that  the  position  y  where  Vo 
vanishes  will  also  be  a  function  of  x.  In  fact  for  numerical  reasons  it  is  more  convenient  to 
treat  <?Ox(x,0)  as  a  third  unknown  to  be  iterated  upon  on  until  qo(x,0)  T  J£°  Tody  ~  0.  The 
numerical  free  boundary  value  problem  specified  by  (3. 7, 3. 8)  can  be  solved  in  principle  by 
adapting  the  procedure  described  in  Hall  and  Lakin  (1988),  however  it  was  found  in  that 
paper  that  the  form  of  the  non-self-similar  solution  can  be  inferred  from  the  self-similar  ones 
by  varying  the  parameter  corresponding  to  G. 

It  remains  for  us  to  discuss  the  nature  of  the  flow  described  above  in  the  neighborhood 
of  the  wall.  We  recall  that  the  solution  which  we  have  obtained  has  u  =  v  =  0,  To  —  T  at 
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the  wall  but  Vo(x,0)  ^  0.  Of  course  the  total  flow  must  satisfy  the  no-slip  conditions  at  the 
wall  so  an  inner  boundary  layer  must  be  present  as  the  wall  is  approached.  We  can  see  from 
(3.3a)  that  in  the  limit  y  — *  0  the  solution  calculeted  above  is  such  that 


u  ~  yu0y(x,0)  -|- - 1-  a~1[[/0(x,Q)E  +  •••]••• 

which  suggests  that  a  new  structure  will  emerge  when  y  falls  to  0(a_1).  Thus  we  define  a 
new  variable  £  by 

C  =  ay 

and  seek  a  solution  for  £  =  0(1).  An  examination  of  the  higher  harmonics  in  (3.3),  which 
are  smaller  than  the  fundamental  for  y  =  0(1),  shows  that  for  £  =  0(1)  all  modes  are 
comparable.  This  is  because  the  cascade  of  energy  from  the  fundamental  down  to  the 
harmonics  is  enhanced  to  such  an  extent  by  the  vertical  diffusion  of  vorticity,  now  comparable 
with  spanwise  diffusion,  that  the  energy  in  the  different  modes  is  of  a  similar  size. 

Thus  in  the  neighborhood  of  the  wall  we  must  replace  (3.3)  by 

u  =  a-1ii(x,  £,  z)  +  0(a~5/3), 
v  =  au(x,£,  z)  +  0(a1/3), 

w  =  aw(x,  £,  z)  -f  0(a-2/3),  (3.9a,b,c,d,e) 

p  =  a2p  +  0(a4/3), 

T  =  T  +  a-1T’(x,  £,  z)  +  0(a~5/3), 


whilst  p  now  expands  as 


P  =  9(*,0  +  °(G  2/3)- 


The  zeroth  order  approximation  to  (2.8)  in  the  wall  layer  then  becomes, 


+  w z 
VU(  -f  wuz 
i>V(  -f  wvz 
VW(  +  wwz 


=  0, 

=  4-  Uzz, 

=  GT  ~  Pt  +  {y<c  +  vzz, 


—  ~Pz  +  +  wZ2, 


vT(  +  wTz 


\  +  f-*} 


k  =  0 


(3.10a, b,c,d,e,f) 
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which  must  be  solved  subject  to 

it  =  v  =  w  —  0, 

u  ->  uOj,(z,0)C  +  [Uo{x,0)E  +  C.C.],  £  -*  oo, 

P  -»  o,  £  -+  oo, 

u  ->  [Vo(x,0)i?  +  C.C.],  C  — ►  oo, 
u>  — *  0,  (  — *  oo, 

q  -*  90(^1°)^  -» oo- 

and  it  is  easy  to  show  that  a  solution  of  the  above  system  can  be  found  by  integrating  from 
£  =  00  to  £  =  0  using  an  appropriate  asymptotic  form  for  £  1.  Thus  the  wall  boundary 

layer  is  passive  even  though  it  is  fully  nonlinear.  Hence  we  can  obtain  the  core  solution  for 
y  =  0(1)  without  reference  to  the  wall  layer  problem  so  we  do  not  discuss  further  the  latter 
problem.  We  shall  now  discuss  self-similar  solutions  of  the  strongly  nonlinear  interaction 
problem  (3.7)  -  (3.8). 

4.  A  Self-Similar  Solution 

We  shall  now  concentrate  on  the  special  case 

ue  =  x1/3,  T  =  x1/3, 

and  introduce  the  similarity  variable  r\  defined  by 

V 

We  then  write 

«o  =  x1/3/'(t/),  v0  =  ^775(7/'  ~  2/},  T0  =  xl/3g{ii),  V'0  =  V(v),  q0  =  x2/3q(r]). 

The  free  boundary  problem  (3.7)  -  (3.8)  can  then  be  written  in  the  simplified  ordinary 


10 


differential  equation  form: 


r+iw-r+i)  = 

V-jU's-  2/s'}  =  -|JW(V'2){ff'V2}' 

H(V2){s~  1  +  x-)  =  0, 

Gcr 

<}'  =  9 

f  =  f  =  0,  T]  =  0,  /'  =  1,9  =  g  =  0,  7/  =  oo. 

In  addition  we  require  that  f,f,q,g  and  g1  are  continuous  at  77  =  rj  where  V  =  0.  The 
problem  can  be  solved  numerically  by  making  guesses  for  1/(0), <j(0),  /"(0)  and  integrating 
(4.1a,b,d)  to  find  V,f,  and  q  with  g  =  1  —  At  the  point  where  V  =  0  we  then  set  the 
terms  proportional  to  V1  in  (4.1a,b)  equal  to  zero  and  integrate  these  equations  together  with 
(4. Id)  to  determine  /, g  and  q  for  77  >  rj.  This  integration  is  carried  out  using  the  computed 
values  of  f,f,f",g,g'  and  q  found  as  q  — »  rj.  For  arbitrary  choices  of  F(0), /"(0),  <7(0)  the 
conditions  /'(oo)  =  1,  5(00)  =  0,  q{ 00)  =  0  will  not  be  satisfied  but  we  can  perform  a 
Newton  iteration  procedure  on  the  wall  values  of  V,f",g'  until  the  conditions  at  infinity  are 
satisfied. 

Numerical  solutions  of  (4.1)  were  in  the  first  instance  obtained  for  a  range  of  values  of  G 
for  S  =  0,  a  =  .2, 1,5.  In  Figure  1  we  show  the  computed  values  of  rj,  K2(0)  and  /"( 0).  We 
see  that  solutions  can  be  obtained  for  G  greater  than  some  finite  value,  in  fact  this  critical 
value  corresponds  to  the  right  hand  branch  of  the  neutral  curve.  At  the  larger  values  of  G 
used  in  the  calculations  the  results  suggest  that  rj  ~  G,  F2(0)  G2,f"(0)  ~  dr-1.  We  shall 

comment  further  on  this  asymptotic  limit  in  the  next  section. 

In  Figures  2,  3,  4  we  show  the  functions  g  and  V2  in  terms  of  77  for  different  values 

of  G  with  a  respectively  equal  to  1,5,  .2.  It  can  be  seen  in  each  case  that  the  boundary  layer 
thickens  as  the  Grashof  number  increases.  In  Figures  2b,  4b  we  note  the  discontinuity  in  /'" 
which  occurs  at  the  transition  layer  position  rj  =  rj.  The  discontinuity  also  occurs  for  the  case 
cr  =  5  shown  in  Figure  3b  but  at  the  values  of  G  used  in  the  calculations  the  discontinuity 
is  not  apparent.  In  each  of  the  calculations  the  function  V2  decreases  monotonically  from 
its  value  at  the  wall  to  zero  at  77  =  rj.  The  temperature  profiles  shown  in  Figures  2c,  3c,  4c 
illustrate  the  large  G  structure  mentioned  above.  Thus  as  G  increases  the  interval  over  which 
g  is  linear  in  77  itself  increases  linearly  with  G.  On  the  other  hand  as  G  is  decreased  towards 
its  linear  critical  value  the  temperature  profile  approaches  it’s  unperturbed  form.  It  is  also 
clear  from  the  calculated  velocity  and  temperature  fields  that  an  important  consequence  of 
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the  mean  flow  being  driven  by  the  vortices  is  that  the  boundary  layer  thickness  is  increased 
from  it’s  unperturbed  value.  More  precisely  we  note  that  in  the  presence  of  vortices  the 
boundary  layer  thickness  is  increased  by  a  factor  G  from  it’s  unperturbed  value. 

A  significant  difference  between  the  calculations  for  the  different  values  of  a  can  be  seen 
in  Figures  2b,  3b,  4b.  We  refer  to  the  fact  that  at  the  smallest  value  of  a,  namely  a  =  .2,  the 
function  f"  has  a  minimum  in  the  region  where  the  vortices  exist  and  a  discontinuity  in  the 
sign  of  its  derivative  across  the  transition  layer.  This  result  is  significant  because  it  means 
that  the  mean  velocity  profiles  associated  with  /'  have  inflection  points  at  the  minima  of 
f"  and  sign  changes  in  the  second  derivative  of  the  mean  downstream  velocity  component 
at  the  transition  layers.  We  stress  that  no  such  points  were  found  for  the  cases  a  =  1,5, 
this  suggests  that  inflection  points  can  only  be  created  by  the  vortices  below  some  critical 
value  of  the  Prandtl  number.  This  means  that  in  low  Prandtl  numbers  flows,  wall  heating 
not  sufficiently  large  to  induce  Taylor-Goldstein  modes  because  5  =  0,  might  still  lead  to 
highly  unstable  inviscid  Rayleigh  waves.  The  secondary  instability  of  the  flows  induced  by 
the  vortices  will  be  discussed  in  the  final  section  of  this  paper. 

In  Figure  5  we  show  /"'  as  a  function  of  tj  for  the  case  a  =  .2.  We  note  that  the  size  of 
the  jump  in  f"  across  the  transition  layer  increases  with  G.  In  addition  when  G  increases  we 
see  that  the  region  between  the  wall  and  the  inflection  point  increases  but  that  remains 
relatively  small  until  the  transition  layer  is  reached.  In  order  to  understand  this  behavior 
it  is  instructive  to  consider  the  limit  G  — >  oo  in  (4.1)  and  see  how  the  underlying  flow 
structure  evolves.  This  limit  will  be  discussed  in  the  next  section,  before  doing  so  we  report 
on  some  calculations  we  have  carried  out  for  the  case  when  buoyancy  effects  are  present.  In 
Figures  6a,b,c  we  show  the  quantities  rjl,  f"(0)  and  V/2(0)  for  the  case  a  =  1,  and  5  =  0, 1,  2. 
The  velocity  and  temperature  profiles  associated  with  the  5  =  2  calculations  are  shown  in 
Figures  6d,e,f,g,  we  note  here  that  the  most  significant  difference  between  the  5  =  0  and 
5  =  2  results  is  that  buoyancy  effects  cause  f  to  overshoot  it’s  free-stream  value  for  a  range 
of  values  of  r/. 

5.  The  Limiting  G  — *  oo  Flow  Structure 

We  stress  here  that,  although  our  asymptotic  analysis  is  carried  out  on  the  similarity 
problem  discussed  in  Section  4,  the  approach  can  be  applied  in  a  similar  way  to  the  nonsimilar 
problem  in  the  limit  x  — ►  oo  with  G  fixed,  see  §7  for  a  brief  discussion  of  this  situation. 

In  order  to  develop  a  large  G  solution  of  (4.1)  it  is  first  convenient  to  write  (4.1a,c)  in 
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the  form 

(t)’  -  £{&-*)r+*fh 

(5.1a,  b) 

-  V)f  + 


which  are  valid  if  V2  >0. 

The  numerical  results  discussed  in  the  previous  section  suggested  that  for  G  1  the 
vortices  are  distributed  over  a  region  of  depth  0(G).  Hence  we  seek  a  solution  of  (5.1)  with 
rj  =  0(Cr).  Actually,  since  G  is  always  multiplied  by  a  when  it  appears,  it  is  more  convenient 
to  treat  crG  as  a  large  parameter  and  define 


r  = 


|{(/'2-l-2//")}  +  ^{(6<7 


l  +  ¥] 


0  = 


Go 


(5.2a) 


It  then  follows  that  the  right  and  left  hand  sides  of  (5.1b)  will  balance  if  /  ~  Gar ,  V  ~  Gcr , 
we  therefore  write 

/  =  G< r#(0)  +  -.., 

(5.2b,  c) 

V  =  GaV{iP)  +  ---, 

and  the  zeroth  order  problem  obtained  from  substituting  the  above  expansions  into  (5.1)  is 


V"  =  (£-)  {cr(ilf/2  —  1  -2®«")  +  #"([1  -0]tf'  +  2«)|?-2, 

\o(T  /  ^  ' 


V2\  1 


(5.3a,  b) 


If  we  assume  that  there  is  no  sublayer  structure  near  0  =  0  then  (5.3)  must  be  solved  subject 
to 

=  0,  V2  =  F2(0),  0  =  0  (5.4) 

and  V(0)  is  a  constant  to  be  determined.  It  can  be  shown  that  for  rj  >  Gcr,  where  no  vortices 
exist,  a  large  (Gcr)  solution  of  (4.1)  can  only  be  developed  if  /'  — >  1,  g  — »  0, 77  — ►  (Gcr)+  so 
in  addition  to  (5.4)  we  require 


V  =  1,  V  =  0,  0  =  1. 


(5.5) 


The  nonlinear  differential  system  specified  by  (5.3),  (5.4),  (5.5)  can  in  general  only  be  solved 
numerically  and  the  solution  will  fix  V'(O)  and  ^"(0).  However  the  expansions  (5.2b,c)  are 
not  uniformly  valid  as  0  — ►  1_  because  in  that  limit  V  ~  (1  —  0)  so  the  constant  term  in 
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the  denominator  of  (5.1b)  will  be  comparable  with  V 2  when  1  —  xp  ~  O(Gcr)-2.  In  order  to 

find  the  appropriate  expansions  in  that  layer  we  observe  that  in  the  neighborhood  of  ip  =  1 

we  can  write  __ 

$  =  xp  +  xp  -f  xp 


V7 2  2  — 

T  =  -3?[1  +  *![^-xl+" 

where  ip  is  a  constant  determined  by  the  solution  in  0  <  xp  <  1  and  xp  is  given  by 


(5.6a,  b) 


xp  =  c(l  —  xp)1+a 


(5.7) 


where  c  is  a  constant.  It  follows  that  in  the  neighborhood  of  xp  =  1  the  expansion  (5.2b) 
must  be  modified  to  give 


/  =  Ga 


0  (0) 
(Ga)2+^ 


(5.8) 


where  <p  =  {Ga)2(  1  -  ip). 


The  function  xp  is  then  found  to  be  given  by 


0=  C 


3 

(b  -f-  - — ~ 

[  Mi  +  VO 


(5.9) 


The  expansion  (5.8)  is  valid  until  0  =  0  where  the  vortex  vanishes.  The  solution  found  above 
is  then  matched  onto  the  solution  of 


r+ j{2/r-/,J  +  i}  =  o 

satisfying  /'( oo)  =  1.  However,  this  outer  layer  is  passive  and  so  we  do  not  pursue  the 
solution  further  here. 

It  follows  from  the  large  Gcr  solution  found  above  that  if  a  <  1  in  this  limit  u"  ~  f"  is 
a  maximum  in  the  layer  near  7]  =  Ga.  In  fact  from  above  it  can  be  shown  that  for  0  =  0(1) 

r  =  -cd + tm* -i)(*+ (c*)2"-”. 

Thus  f"  becomes  large  if  a  <  1,  this  is  consistent  with  our  numerical  work.  Numerical 
solutions  of  (5.3)  suggest  that  the  constant  c  is  always  positive  so  that  f"  is  positive  or 
negative  near  ?/  =  Ga  depending  on  whether  or  not  a  is  less  or  greater  than  1.  Moreover,  in 
the  case  a  <  1,  f"  attains  its  largest  value  near  77  =  Ga  and  this  value  ~  (Gcr)2(1“<7h  Since 
we  can  show  that  the  solution  for  7 ]  >  G0  must  have  f'"  <  0  we  see  that  the  results  shown 
in  Figure  4  which  showed  a  jump  in  the  sign  of  f"  at  rj  are  confirmed  and  indeed  a  similar 
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result  would  be  found  for  any  value  of  a  <  1.  Furthermore  since  a  small  rj  solution  of  (4.1) 
shows  that  /"'  is  positive  near  r;  =  0,  the  continuity  of  /"'  in  (0,7;)  implies  the  existence  of 
an  inflection  point  in  the  velocity  profile  for  any  Prandtl  number  less  than  unity. 

Some  numerical  solutions  of  (5.3)  are  shown  in  Figures  7,  8,  9.  Figure  7  shows  ^"(0) 
and  V(0)  as  a  function  of  a.  The  computed  velocity  field  in  the  large  Ga  limit  similarly 
agrees  well  with  the  full  solutions.  Of  particular  interest  is  the  dependence  of  the  function 
'i!'"  on  ip  and  a.  In  Figures  8a, b  we  show  plots  of  and  fl/'"  as  functions  of  ip  and  a.  The 
results  shown  confirm  that  inflectional  velocity  profiles  exist  only  for  a  <  1.  At  finite  values 
of  Go  our  calculations  suggest  again  that  inflectional  profiles  exist  only  for  o  <  1  though  an 
exhaustive  check  of  this  possibility  in  the  a  —  G  plane  was  not  carried  out.  In  the  case  <7^1 
the  results  of  Figure  7  suggest  that  for  any  a  >  1,  'I'"'  is  always  negative  whilst  for  o  <  1  ty'" 
is  positive  in  an  interval  (ip’,l)  for  some  ipm  <  1.  In  order  to  demonstrate  the  behavior  of 
W  near  o  =  1  we  have  in  Figure  9  plotted  results  for  the  cases  a  —  .98,  .99, 1, 1.01, 1.02.  We 
see  that  tends  to  a  negative  constant  as  'I'  — ♦  1  for  o  =  1,  whilst  for  o  slightly  less  than 
1  W  becomes  large  and  positive  in  a  small  interval  near  ip  —  1.  For  o  slightly  greater  than 
1  the  function  'I'"'  is  always  negative  but  increases  significantly  in  magnitude  near  ip  =  1. 
We  note  that  the  finite  value  of  for  o  =  1  occurs  because  the  next  correction  term  in 
(5.6a)  when  a  =  1  is  proportional  to  ip3.  The  implications  of  the  above  results  for  secondary 
instability  theory  will  be  discussed  in  §6. 


6.  Rayleigh  or  Taylor-Goldstein  Inviscid  Breakdown  Induced  by  Streamwise 
Vortices 


Now  let  us  consider  the  inviscid  instability  of  the  finite  amplitude  vortex  structures 
described  in  the  previous  section.  The  approach  we  take  is  based  on  the  work  of  Hall  and 
Horseman  (1991)  on  the  instability  of  Gortler  vortices.  The  most  important  property  of 
an  inviscid  instability  in  a  boundary  layer  flow  is  that  it  operates  on  the  same  streamwise 
lengthscale  as  the  boundary  layer  thickness.  In  terms  of  (2.8)  this  means  that  for  the  inviscid 
wave  disturbance  ^  ~  R}^2.  In  addition  inviscid  waves  are  time-dependent  so  that  terms 
lit,  Tt  must  be  added  to  (2.8b,c).  Thus  we  now  write 


{  u  ^ 

'  u{x,y,z)  > 

(  U{x,y,z)  \ 

V 

v{x,y,z) 

R1/2V(x,y,z) 

IV 

= 

W{x,y,  z) 

!t'!HV{x,yyz) 

p 

p(x,y,z) 

RP{x ,  y,  z) 

\T  ) 

\  T(x,y,z )  ) 

\  0{x,y,z)  / 

(6.1) 


Here  the  first  term  corresponds  to  the  combined  mean  flow-vortex  state  driven  by  buoyancy 
effects  whilst  the  second  term  represents  an  inviscid  travelling  wwe  disturbance  of  arbitrarily 
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small  amplitude  A.  We  again  denote  the  buoyancy  parameter  S  =  G/R  and  then  substi¬ 
tute  (6.1)  into  (2.8)  and  equate  terms  of  order  A.  In  the  limit  R  — >  oo  the  zeroth  order 
approximation  to  this  system  is 


iaV  +  Vy  +  W,  =  0, 

ia(u  -  c)U  +  Vuy  -f-  Wuz  =  —iaP, 


ia(u  —  c)V  =  —  Py  ■+  Sd, 


ia(u  —  c)W  =  —Pz, 
ia{u-c)6  +  VTy  +  WTz  =  0. 


(6.2a,b,c,d,e) 


Here  we  have  replaced  ^  by  the  wavespeed  c  and  the  appropriate  boundary  conditions  are 
V  =  0,  y  =  0,  oo.  With  the  function  u  specified  by  the  steady  nonlinear  vortex  problem 
discussed  earlier  the  system  (6.2)  and  the  boundary  conditions  constitute  an  eigenvalue 
problem 

a  =  a(c ,  S). 


Note  that  if  the  basic  state  is  dependent  only  on  x,y  then  (6.2)  can  be  reduced  to 


qt'  _ 

(S  -c)(V„-a‘‘V)-u„V  =  -—V,  (6.3) 

which  is  the  so-called  Taylor- Goldstein  equation.  If  buoyancy  forces  are  negligible,  5  =  0, 
(6.3)  reduces  to  Rayleigh’s  equation.  However  when  u  is  a  function  of  z  no  simple  general¬ 
ization  of  (6.3)  is  available,  but  we  can  eliminate  U,  V,  W  and  9  from  (6.2)  to  give 


(  P ,  \  °2p  c.  [  P,Tt  +  ft T, 

\Ju  —  c]2  /  v  \Ju  —  c\2  J  z  [u  —  c]2  [u  —  c]2[5Ty  —  q2(u  —  c)2] 


(6.4) 


If  5  =  0  the  above  equation  reduces  to  the  equation  obtained  by  Hall  and  Horseman  (1991) 
in  their  discussion  of  inviscid  instabilities  induced  by  Gortler  vortices.  The  solution  of  (6.4) 
is  of  course  a  nontrivial  task  bearing  in  mind  the  fact  that  u  must  in  general  be  determined 
numerically.  Here  we  shall  concentrate  on  the  case  when  5  =  0  so  that  we  are  in  effect 
limiting  our  analysis  to  the  determination  of  whether  vortex  instabilities  induced  by  wall 
heating  can  trigger  a  rapidly  growing  Rayleigh  instability.  Thus  we  shall  now  confine  our 
attention  to  the  solution  of  the  eigenvalue  problem 

Pz 

[u  -  c]2 

Py  =  0,  y  =  0,  P  -*•  0,  y  —>  oc. 


cPP 
[u  —  Cp 


=  0 


(6.5) 
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with  ti(x,  y,z)  determined  by  the  nonlinear  vortex  equations  in  the  presence  of  wall  heating 
but  with  5  =  0.  In  particular  we  will  investigate  the  case  when  u  corresponds  to  the  small 
wavelength  solution  discussed  in  §3,4.  We  recall  that  in  that  limit  vortices  are  confined  to  a 
region  of  depth  0(1)  adjacent  to  the  wall  with  boundary  layers  of  thickness  a-1,  a-2/3  at  the 
wall  and  at  the  edge  of  the  region  of  vortex  activity  respectively.  Above  the  a~ 2/3  transition 
layer  the  flow  is  determined  by  the  unperturbed  boundary  layer  equations  and  the  mean 
parts  of  u,  uy  with  respect  to  z  are  continuous  across  the  layer.  It  follows  that  a  solution 
of  (6.5)  can  be  sought  with  no  z  dependence  where  vortices  are  absent.  In  the  lower  part  of 
the  boundary  layer,  where  u  expands  as  in  (3.3a),  we  find  that  P  takes  the  form 

P  =  Pq  +  a-1  Pi  +  a~2P2  -f  a~^P3  cos  az  +  a~3P4  +  •  •  •  (6.6) 

where  Pa,  Pi  etc.  depend  only  on  x,y.  If  we  expand  a,c  in  the  form 

(a,c)  =  (q0,co)  +  a"1(a1,c1)  +  (6.7) 


then  we  find  that  Pq,  P3  satisfy 


Poy  \  ft2P0 
[uo-CoP/j,  [uo-co]2 


(6.8) 


and 

P3  =  -2(fi0  -  Co)2  (U0P0y[uo  ~  c0]-3)y  •  (6.9) 

Equation  (6.8)  is  simply  the  Rayleigh  equation  for  a  unidirectional  flow  u0  so  that  the  vortex 
does  not  have  a  direct  effect  on  the  zeroth  order  Rayleigh  problem.  However  it  does  have  a 
significant  indirect  effect  because  uo  is  determined  by  the  vortices.  In  the  region  above  the 
a-2/3  transition  layer  (6.8)  again  gives  the  correct  zeroth  order  approximation  to  the  inviscid 
stability  problem.  Across  the  layer  xio,  uoy  are  continuous  whilst  uoyy  is  discontinuous. 
In  fact  this  discontinuity  in  u0yy  is  smoothed  out  within  the  transition  layer  by  viscous 
effects.  An  examination  of  (6.5)  in  the  transition  layer  shows  that  as  long  as  c0  /  u0  in 
this  layer  then  Po,  Poy  are  continuous  across  the  layer.  Similarly  the  a-1  wall  layer  is  passive 
so  that  if  we  are  not  concerned  with  neutral  waves  propagating  downstream  with  the  mean 
downstream  fluid  speed  in  the  transition  layer  then  it  is  sufficient  fur  us  to  solve  (6.8)  with  u0 
determined  by  the  vortices  in  0  <  y  <  y  and  by  the  boundary  layer  equations  for  y  >  y.  For 
definiteness  we  consider  the  case  when  u0  has  the  similarity  form  discussed  in  §4.  Thus  we 
write  u0  =  xl^3f'(r})  with  /  determined  by  (4.1),  it  is  convenient  to  rescale  Qo,Co  by  writing 


q0  =  cr0x~l/3,  cq  =  c0xl/3 
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so  that  the  eigenvalue  problem  d0  =  c»o(coi  G )  becomes 

/  Pon  \  _  &'jPo  _  ,, 

Vi 

Pqt)  =  0,  rj  =  0,  Pq  — >  0,  ij  —>  x.  (6.10) 

In  fact  the  similarity  solution  we  have  considered  is  particularly  important  in  vortex-wave  in¬ 
teraction  theory.  In  the  latter  theory,  see  Hall  and  Smith  (1991),  a  small  amplitude  Rayleigh 
or  Tollmien-Schlichting  wave  system  interacts  with  itself  to  drive  a  strong  vortex  field  which 
itself  acts  back  on  the  wavefield.  In  the  present  context  a  Rayleigh-vortex  interaction  occurs 
if  A  in  (6.1)  is  choosen  appropriately.  In  this  case  the  wave  system  drives  the  vortex  in  the 
critical  layer  and  the  fact  that  the  wave  system  must  remain  neutral  as  it  moves  downstream 
means  that  similarity  solutions  of  the  inviscid  equation  describing  the  wave  are  possible  only 
if  n  =  |. 

In  the  absence  of  wall  heating  the  basic  state  has  f  determined  by  the  Falkner-Skan 
problems  for  a  pressure  gradient  proportional  to  a.--1/3,  and  since  the  velocity  profile  is 
non-inflectional  no  unstable  inviscid  eigenvalues  exist.  We  saw  in  the  previous  section  that 
inflectional  streamwise  velocity  profiles  are  generated  by  wall  heating  whenever  the  Prandtl 
number  ct  is  less  than  unity. 

In  Figures  10,  11,  12,  13  we  show  results  we  have  obtained  by  solving  (6.10)  for  the  case 
a  —  .2  and  a  range  of  values  of  G.  In  our  calculations  we  have  kept  do  real  and  computed 
the  corresponding  complex  value  of  co-  We  see  in  Figure  10  that  at  each  value  of  G  there  is 
a  band  of  unstable  modes  present  to  the  right  of  a  finite  value  of  d0  .  Each  mode  becomes 
neutral  and  disappears  when  the  wavespeed  is  equal  to  the  fluid  speed  at  the  inflection  point 
of  the  velocity  profile.  Note  here  that  all  our  calculations  were  for  cases  where  such  an 
inflection  point  exists,  at  sufficiently  small  values  of  G  the  basic  mean  profile  approaches 
its  unperturbed  values  and  no  inflection  points  exist.  The  unstable  mode  persists  for  all  do 
greater  than  the  neutral  value  but  has  (doCo)t  tending  to  zero  as  do  —*  oo,  in  this  limit  cor 
approaches  the  fluid  speed  at  the  discontinuity  in  /"'  at  the  transition  layer.  An  analysis 
of  that  limit  shows  that  the  present  analysis  breaks  down  when  do  =  0(loga).  In  this 
wavenumber  regime  the  Rayleigh  wave  has  a  two  layer  structure  (of  depth  a-2^3,  (logo)-1  ) 
at  the  transition  layer  and  is  effectively  zero  elsewrhere.  We  do  not  give  details  of  the  behavior 
for  do  ~  |  logn|  because  the  mode  is  neutral  there  and  the  most  unstable  growth  rates  occur 
for  do  =  0(1).  We  see  that  the  wavenumber  of  the  most  unstable  mode  increases  as  the 
Grashof  number  increases,  this  is  of  course  due  to  the  thickening  of  the  boundary  layer. 
Some  eigenfunctions  associated  with  the  modes  are  shown  in  Figures  12,  13.  We  see  the 
concentration  near  rj  of  the  eigenfunction  of  the  second  mode  as  its  growth  rate  approaches 
zero  at  large  do. 


7.  Conclusion 


We  have  shown  that  wall  heating  produces  large  amplitude  streamwise  vortex  structures 
which  completely  alter  the  boundary  layer  in  which  they  develop.  We  have  concentrated  our 
attention  on  self-similar  flows  which  enabled  us  to  solve  the  mean  flow-vortex  interaction 
problem  by  reducing  it  to  a  set  of  nonlinear  differential  equations.  A  significant  result  which 
we  found  was  that  the  mean-state  modified  at  zeroth  order  can  have  inflection  points  whereas 
the  unperturbed  state  does  not.  For  the  special  case  when  the  driving  pressure  gradient  is 
proportional  to  we  found  that  inflectional  profiles  are  only  generated  when  the  Prandtl 
number  is  less  than  unity. 

The  importance  of  the  inflectional  profiles  is  that  they  are  highly  unstable  to  Rayleigh 
waves  growing  on  a  streamwise  lengthscale  0 shorter  than  that  over  which  the  near 
state  develops.  This  means  that  a  boundary  layer  inviscidly  stable  in  the  absence  of  wall 
heating  can  be  made  massively  inviscidly  unstable  by  streamwise  vortex  restructuring  of 
the  boundary  layer.  It  is  of  course  relevant  to  question  whether  the  significantly  different 
mean  flow  character  corresponding  to  the  cases  a  <  1,  cr  >  1  is  a  function  of  the  particular 
similarity  flow  considered.  In  order  answer  this  question  we  note  that  the  more  general  form 
of  the  similarity  solution  given  by  (4.1)  with  S  =  0  is 

n  — 1 

v  =  yx  2  , 


uo  =  xnf'{r]),v  o  = 


Pn-1 


rn-l_„  .  n  +  l  r 
—n!  +  — / 


Vo  =  V(v),T0  =  x'^giri), 

f“  +  {  —-SI"  -  nf 2  +  n)  =  -I{/V,}'//(V'a), 

V  + 


(7.1) 


/  =  /'  =  0,7  =  0,/'  =  1,5  =  0,7?  =  oc, 

/,  /',  5,  ^continuous  at  7?  =  Tj,  where  V' (77)  =  0. 

In  order  to  see  whether  inflectional  profiles  exist  at  large  Go  we  can  repeat  the  analysis  of  §5 
by  seeking  a  solution  for  0.  The  expansion  procedure  follows  exactly  that  as  of  §5  and  the 
key  result  is  that  the  correction  term  ip  in  (5.6a)  is  independent  of  n  and  again  proportional 
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to  (— ip  +  l)£r+l.  This  means  that  /"'  becomes  large  as  tp  — >  1_  and  changes  sign  when  a 
passes  through  1.  This  result  can  be  used  to  infer  that  /"'  is  positive  as  xjj  — >  1—  for  a  <  1 
and  coupled  with  the  fact  that  /"'  must  be  negative  for  //  >  Tj  and  for  //  <C  1  we  find  that  for 
any  n  the  large  Go  limit  leads  to  inflectional  profiles  for  a  <  1.  This  argument  suggests  that 
the  results  obtained  for  the  special  case  n  =  1/3  are  typical  and  that  inflectional  profiles  are 
produced  whenever  o  <  1  at  sufficiently  large  values  of  Go.  In  fact  the  large  Grashof  number 
analysis  of  §5  can  be  reformulated  as  a  large  x  asymptotic  solution  of  the  full  interactive 
problem.  This  can  be  done  for  ue  ~  xn,T  ~  xm  for  any  positive  m  and  it  is  found  that 
inflectional  profiles  are  again  only  created  for  Prandtl  numbers  less  than  unity.  This  suggests 
that  all  heated  boundary  layer  flows  of  fluids  with  Prandtl  number  less  than  unity  are  caused 
to  become  inviscidiy  unstable  by  a  streamwise  vortex  restructuring  of  the  flow. 
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Figure  lc.  The  quantities  /"( 0),  \/2(0),  rj  for  the  similarity  solution  with  a  =  .2,  5  =  0. 
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Figure  2a.  The  function  /'(j?)  for  S  =  0,  <7=1,  G  =  2.,  4.,  6.,  8. 
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Figure  3a.  The  function  /'(t?)  for  S  =  0,  <7  =  5,  G  =  .3, 1.1, 1.9, 2.7. 
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Figure  6d.  The  function  /'(tj)  for  a  =  1,  5  =  2  and  G  =  1.5, 3, 4. 5, 6. 
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Figure  13.  The  imaginary  part  of  the  eigenfunction  V0  =  as  a  function  of  rj  for 

5  =  0,  a  =  .2,  G  =  28.83  and  d0  =  .1,  .3,  .5,  .7,  .9, 1.1, 1.3, 1.5, 1.7, 1.9. 
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